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Nomenclature 

a = 
D = 
E \> £3 = 


n = 

Q = 
T = 



<f> = 
* = 

Subscripts 

/, o = 
r = 
1,2 = 


absorption coefficient of layer, m~ 
thickness of plane layer, m 
exponential integral functions, 

E m (*)= ( ii m - 2 exp(-x/fi)dn 

J o 

index of refraction 
heat flux, W/m 2 
absolute temperature, K 
temperatures of surrounding radiating re- 
gions, K 

coordinate normal to boundary of plane 
layer, m; X-x/D 

optical depth, (a + o s )x\ k d , optical thickness, 
(a + o s )D 

Stefan-Boltzmann constant, W/m 2 -K 4 
scattering coefficient of layer, m -1 
reflectivity of interface for internally incident 
radiation 

transmissivity of surface for externally inci- 
dent radiation 

dimensionless temperature distribution 
dimensionless radiative heat flux 


incoming and outgoing radiation 
radiative quantity 

the hotter and cooler surroundings of the 
layer 


Introduction 

The index of refraction can considerably influence the tem- 
perature distribution and radiative heat flow in semitranspar- 
ent materials such as some ceramics. For external radiant 
heating, the refractive index influences the amount of energy 
transmitted into the interior of the material. Emission within 
a material depends on the square of its refractive index, and 
hence this emission can be many times that for a blackbody 
radiating into a vacuum. Since radiation exiting through an 
interface into a vacuum cannot exceed that of a blackbody, 
there is extensive reflection at the internal surface of an in- 
terface, mostly by total internal reflection. This redistributes 
energy within the layer and tends to make its temperature 
distribution more uniform. 

The purpose of the present analysis is to show that, for 
radiative equilibrium in a gray layer with diffuse interfaces, 
the temperature distribution and radiative heat flux for any 
index of refraction can be obtained very simply from the results 
for an index of refraction of unity. For the situation studied 
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here, the layer is subjected to external radiative heating incident 
on each of its surfaces. The material emits, absorbs, and iso- 
tropically scatters radiation. For simplicity the index of re- 
fraction is unity in the medium surrounding the layer. 

The surfaces of the layer are assumed diffuse. This is prob : 
ably a reasonable approximation for a ceramic layer that has 
not been polished. When transmitted radiation or radiation 
emitted from the interior reaches the inner surface of an in- 
terface, the radiation is diffused and some of it thereby placed 
into angular directions for which there is total internal reflec- 
tion. This provides a trapping effect for retaining energy within 
the layer and tends to equalize its temperature distribution. 

An analysis of temperature distributions in absorbing-emit- 
ting layers, including index of refraction effects, was developed 
by Gardon (1958) to predict cooling and heat treating of glass 
plates. The interfaces were optically smooth; the resulting spec- 
ular reflections were computed from the Fresnel reflection 
laws. This provides a somewhat different behavior than for 
diffuse interfaces. A similar application was for heating that 
occurs in a window of a re-entry vehicle (Fowle et al., 1969). 
A number of recent papers (Rokhsaz and Dougherty, 1989; 
Ping and Lallemand, 1989; Crosbie and Shieh, 1990) further 
examined the effects of Fresnel boundary reflections and non- 
unity refractive index. Other examples of analyses of both 
steady and transient heat transfer to single or multiple plane 
layers (Amlin and Korpela, 1979; Tarshis et al., 1969) have 
used diffuse assumptions at the interfaces as in the present 
study. 

Analysis 

A plane layer of ceramic material has thickness D as shown 
in Fig. 1 . It absorbs, emits, and isotropically scatters radiation. 
The limiting case is considered here where the layer temperature 
distribution is dominated by radiation so heat conduction is 
neglected. The material has a constant index of refraction; it 
is the effect of the index of refraction that is investigated here. 
Because of the large amount of scattering in many ceramic 
materials, the interfaces between the ceramic and the sur- 
rounding air or vacuum are assumed to be diffuse. The re- 
fractive index of the surroundings is unity. As shown in Fig. 
1, the layer is subjected to radiation from the surroundings at 
T sl and T sl on the two sides x=0 and x = D. It is assumed that 
the surroundings act as black environments so the incident 
energies on the two sides are <? H (0) = aF 4 ! and q r i(D) =oT * 2 . 
Inside the layer there are outgoing and incoming fluxes, q 0 
and q h at each interior surface, as shown on the figure. Since 
scattering is included, the local optical depth is related to the 
x coordinate by k = (a + a 5 )x. 

The temperature distribution inside the layer is governed by 
the integral equation given by Siegel and Howell (1981), mod- 
ified with the index of refraction factor as 



Fig. 1 Layer geometry, coordinate system, and nomenclature of heat 
fluxes at boundaries 
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( 10 ) 


rt : orV) = ^ (<7„(0)£ 2 (k) +^„(k d )£;(k- 0 -k)] 
n 2 r° 

+ — ar 4 (K*)£,(U-K*l)t/ K * (1) 

2 J 0 

The radiative flux, which is a constant through the layer for 
the present conditions of radiative equilibrium, can be found 
from the temperature distribution as (Siegel and Howell, 1981), 

Q r = q o (0)-2q o (K D )E 3 (K D )-2n 2 \ D oT A (K)E 2 (K)dK (2) 

J o 

The following dimensionless groups are now defined: 


/» 2 g7- 4 (K)-<7 0 (>c 0 ) 

<7 o (0) - q 0 (K D ) 


(3a) 


¥ = - 


m 


q 0 ( 0) - <7o(kd) 

Equations (1) and (2) then become (note that X= x/D = k/k d ) 


$(X')E\(k d \X- X* I )dX* (4) 

2 2 J 0 


¥=1-2*4 ^(X)E 2 (K D X)dX 


(5) 


In Eqs. (4) and (5) ¥ and ¥ are not functions of n. Hence 
to obtain ¥ and ¥ for all n> 1 it is .necessary to solve Eq. (4) 
only once for each k d and use each result to evaluate Eq. (5). 
For the special case when the index of refraction of the layer 
is 1 , there are no reflections at the interfaces so that q o (0) = oT A u 
Qo(kd) ~ °T a s2 , and the dimensionless groups become 

T a (k)-T a s2 

4>W = — =4 — rj— ; (6a) 

* sl “ * s2 


¥ = 


art -art 


m 


For n > 1 , however, the ¥ and ¥ in Eq. (3) contain the outgoing 
boundary fluxes q Q ( 0) and q Q (K D ) that are unknown, so the 
solution has not yet provided T(k) and q r for the cases when 
n > 1 . In order to find these quantities the <7 o (0) and q 0 {K D ) 
must be obtained in terms of known quantities. This is ac- 
complished by looking at the boundary conditions in detail. 

At the diffuse interfaces the internal fluxes are related to 
the transmission of external flux and the reflection of internal 
flux by 

<7o(0) = oT a s \t° + <?,(0)p' (7a) 

q 0 (K D ) = aTt 2 T° + q l {K D )p i (lb) 


At the inside surfaces of the two boundaries there are the 
following relations between the radiative flux and the outgoing 
and incoming fluxes: 


Qr(0) = q o (0)-q,{0) (8 a) 

qr(*D)= - Qo(*d) + Q&d) ( 86 ) 

Equations (la) and (8tf) are combined to eliminate <7,{0) and 
similarly for Eqs. (76) and (86) to eliminate qfa D ). This yields 


<7»(0) = t~ ;■ (t°°Tsi ~ Qrf>‘) 

1 “P 

( 9a) 

7„(k d ) = --U (t°oT a 2 + q^) 

1 -p 

( 9b) 


These relations are substituted into Eq. (36) to eliminate the 
q Q ' s. This yields the radiative flux for any index of refraction 
in terms of the value of ¥. Use is also made of the relation at 
an interface (Richmond, 1963) that r° = (l - p)n 2 . The result 
for q f Jor any n is, 


Qr 

o(T A st -T A s2 ) 


rr¥ 




i-p' 


Following the same procedure, the temperature distribution is 
found by starting with Eq. (3ff), using Eq. (9) to eliminate the 
q Q 's and then using Eq. (10) to eliminate q r . This yields the 
T(X) for any n as 


T*(X) - T& 
71, -T* a 


*(X) + 


1-p' 




i-p' 


(ii) 


To use these relations values of p are needed for various 
refractive indexes. The externally incident radiation is diffuse. 
Although the interfaces are probably not optically smooth, it 
is assumed that each bit of roughness acts as a smooth facet 
so that the reflectivity can be obtained from the usually used 
interface relations for a nonabsorbing dielectric medium. Then 
by integrating the reflected energy over all incident directions 
the relation for p\ri) is (Richmond, 1963), 


. . ifl (3«+lX»-D n 2 (n 2 - 1) 2 (n-\\ 

<"»- 1 -p(i — ~<? 7 Tr w 

2w 3 (w 2 + 2p— 1) 8« 4 (« 4 +l) . . ■) 


As discussed by Cox (1965), in the fairly transparent spectral 
regions of ceramic materials, the extinction coefficient in the 
complex index of refraction is usually not large enough to 
affect the surface reflectivity significantly, so Eq. (12) for 
nonattenuating dielectrics gives good results. The extinction 


Table 1 Dimensionless temperature distribution, $ 


X 

Optical Thickness, k d jj 

0.1 

0.3 

1.0 

3.0 

10 

30 

100 

0 

.5710 

.6419 

.7582 

.8693 

.9495 

.9819 

.9948 

.05 

.5619 

.6229 

.7230 

.8211 

.8966 

.9304 

.9456 

.10 

.5541 

.6072 

.6946 

.7819 

.8511 

.8828 

.8974 

.20 

.5397 

.5786 

.6429 

.7088 

.7627 

.7876 

.7994 j 

.30 

.5262 

.5517 

.5942 

.6384 

.6750 

.6920 

.6999 

.40 

.5130 

.5257 

.5468 

.5690 

.5874 

.5961 

.6002 

.50 

.5000 

.5000 

.5000 

.5000 

.5000 

.5000 

.5000 

.60 

.4870 

.4743 

.4532 

.4310 

.4125 

.4038 

.4002 

.70 

.4738 

.4483 

.4058 

.3616 

.3249 

.3079 

.3004 

.80 

.4603 

.4214 

.3571 

.2912 

.2372 

.2123 

.2007 

.90 

.4459 

.3928 

.3054 

.2181 

.1488 

.1171 

.1026 

.95 

.4382 

.3771 

.2770 

.1789 

.1034 

.0695 

.0543 

1.0 

.4290 

.3581 

.2419 

.1307 

.0505 

.0181 

.0052 


Table 2 Dimensionless heat flux, if 


Optical 

Thickness 


0.1 

.9157 

0.3 

.7936 

1.0 

.5534 

3.0 

.3017 

10 

.1168 

30 

.0419 

100 

.0122 
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Dimensionless position, X = x/D 
(a) 



Dimensionless position, X = x/D 
ib) 



Dimensionless position. X = x/D 
(c) 

Fig. 2 Effect of Index of refraction on dimensionless temperature dis- 
tributions for various optical thicknesses: (a) index of refraction, n = 1; 
(b) Index of refraction, n = 2; (c) index of refraction, n = 4 



Optical thickness, k 0 = (a + < 7 S )D 

Fig. 3 Effect of index of refraction on radiative heat flux through the 
layer as a function of optical thickness 


effect can be important in highly absorbing wavelength bands. 
General property relations as given by Hering and Smith (1968) 
would then be required to estimate average surface properties 
for use in the present gray analysis. 

Equations (10)— (12) provide the means for obtaining the 
temperature distributions and radiative heat fluxes for a layer 
with n> 1 from the results for n = 1. The solutions of Eqs. (4) 
and (5) for $ and * were obtained in another context by Heaslet 
and Warming (1965). In connection with some of our work 
(Siegel, 1987), a computer program was available; solutions 
were obtained from that program to have results in a conven- 
ient tabular form for the desired optical thicknesses. As dis- 
cussed in that work the program uses a Gaussian integration 
subroutine, and 80 grid points were used across the layer. Cubic 
spline interpolation is used to obtain the irregularly spaced 
integration points for the Gaussian routine. To deal with the 
singularity £,( 0)-oo, the integral in Eq. (4) was evaluated 
analytically for X* very close to /Tby assuming 4>(A"*) ~${X) 
so that $(X) could be taken out of the integral for very small 
\X-X* I (see Siegel, 1987). The results were in agreement with 
those of Heaslet and Warming (1965). 


Results and Discussion 

The and * were evaluated from Eqs. (4) and (5). The 
results are in excellent agreement with the plotted values of 
Heaslet and Warming (1965) and are given for convenience in 
Tables 1 and 2 for 0.1 <k d < 100. Equations (10) and (11) were 
then used to obtain dimensionless radiative heat fluxes and 
temperatures for n> 1. The required p 1 were obtained from 
Eq. (12). 

Figure 2 shows the dimensionless temperature distributions 
for n = 1 , 2, and 4. For small optical thickness the dimensionless 
distribution approaches 0.5, and as k d —oo the profile becomes 
linear extending from 1.0 to 0. The effect of increasing n is 
to decrease the range of the temperature distributions, and for 
a fixed k d they move closer to 0.5 as n is made larger. The 
dimensionless profiles are all rather linear. The fact that the 
profiles become more uniform with increasing n is the result 
of increasing internal reflections within the absorbing and scat- 
tering layer and decreasing transmission into the layer. Since 
each element in the layer is in radiative equilibrium all locally 
absorbed radiation must be reemitted. Since scattering is as- 
sumed isotropic, local scattering is added to the local emission. 
The large amount of internal reflection tends to equalize the 
energy throughout the layer and flatten the temperature dis- 
tributions. The effect of n is calculated very easily from the 
simple relation given by Eq. (11), where p l is a function of n. 

The effect of n and k d on the dimensionless radiative heat 
flux through the layer is in Fig. 3. The heat flux decreases as 
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the optical thickness increases, and the effect of n is quite 
pronounced in altering the heat flux. The curve becomes very 
flat for n = 4. This is because interface reflections and increased 
internal emission (as a result of the n 2 factor) have become 
quite strong in regulating the heat transfer. The effect of optical 
thickness is thereby suppressed. At large optical thicknesses 
this increases the radiative flux as compared with the results 
for 1. As the layer becomes transparent k d — 0, ^ — 1 and 
the dimensionless flux from Eq. (10) approaches n 2 (\ - p 1 )'/ 
0+p'). 
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The Prandtl Number Effect on Melting Dominated by 
Natural Convection 


J. S. Lim 1 and A. Bejan 1 * 2 


Nomenclature 

B = constant, Eq. (13) 
c P = liquid specific heat 
Fo = Fourier number = at/H 2 
g = gravitational acceleration 
Gr = Grashof number = g(3ATH*/p 2 
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h sf = latent heat of melting 
H - height 

k = liquid thermal conductivity 
Nu = Nusselt number = q' (kAT) 

Pr = Prandtl number = p/a 

q' - overall heat transfer rate, W/m 

Ra = Rayleigh number = g&ATH 2 /(ap) 

Ra ; = Rayleigh number = g(3ATz 3 /(ap) 
s = thickness of conduction zone 
Ste = Stefan number — c P AT/h s j 
T h = temperature of isothermally heated wall 
T m = melting point 

AT = temperature difference -T h —T m 
z = height of convection zone 
a = liquid thermal diffusivity 

0 = liquid coefficient of volumetric thermal expansion 
6 = thermal boundary layer thickness 
v = kinematic viscosity 
() c = critical 


Objective 

In a recent paper, Gobin and Benard (1990) considered the 
task of correlating the heat transfer data for melting in the 
presence of natural convection when the Pr value of the liquid 
phase is considerably smaller than 1. Earlier correlations were 
reported by Webb and Viskanta (1986), and Beckermann and 
Viskanta (1989). Additional low-Pr studies were conducted by 
Webb and Viskanta (1985), Gau and Viskanta (1986), Wolff 
and Viskanta (1987), and Beckermann (1989). To correlate the 
low-Pr data is an important and timely task, especially in view 
of the voluminous work that has been dedicated to situations 
in which Pr is greater than 1 . 

For the convection-dominated regime known also as quasi- 
stationary melting, Gobin and Benard (1990) correlated their 
low-Pr numerical data with the formula: 

Nu = 0.29Ra°' 27 Pr 018 (1) 

They noted that this correlation does not agree with the 
Nu~(RaPr) 1/4 trend that might be expected from the single- 
phase natural convection scales for low Prandtl numbers (Be- 
jan, 1984). They concluded that: 

1 The relevance of the group (RaPr) is not verified by their 
numerical results for convection-dominated melting, and 

2 Further work is required to determine the scaling laws 
that govern the transition from the initial (conduction) regime 
to the final (convection) regime of the process of melting by 
side heating. 

These two conclusions defined the work presented in this 
note. In it we report the correct scales of natural convection 
melting when the Prandtl number is small. We then construct 
a scaling-correct heat transfer correlation that spans the entire 
range of Prandtl numbers. 

Scale Analysis 

The scales of the natural convection melting process can be 
determined by extending Jany and Bejan’s (1988) high-Pr the- 
ory to the range of low Prandtl numbers represented by liquid 
metals. The theory is based on the geometric fact that during 
the transition from conduction to convection-dominated melt- 
ing the melt region contains two distinct zones. As shown in 
Fig. 1 , the upper zone of height c is ruled by convection (namely, 
distinct boundary layers, 6), while the lower zone of uniform 
thickness s and height ( H-z ) is governed by horizontal con- 
duction. It is assumed that the flow is laminar, and that 
Ste« 1. 

The conduction thickness is well known, 

s~H( SteFo) l/2 (2) 
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